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We theoretically investigate the momentum-resolved radio-frequency spectroscopy of a harmon- 
ically trapped atomic Fermi gas near a Feshbach resonance in the presence of equal Rashba and 
Dresselhaus spin-orbit coupling. The system is qualitatively modeled as an ideal gas mixture of 
atoms and molecules, in which the properties of molecules, such as the wavefunction, binding en- 
ergy and effective mass, are determined from the two-particle solution of two-interacting atoms. 
We calculate separately the radio-frequency response from atoms and molecules at finite tempera- 
tures by using the standard Fermi golden rule, and take into account the effect of harmonic traps 
within local density approximation. The total radio-frequency spectroscopy is discussed, as func- 
tions of temperature and spin-orbit coupling strength. Our results give a qualitative picture of 
radio-frequency spectroscopy of a resonantly interacting spin-orbit coupled Fermi gas and can be 
directly tested in atomic Fermi gases of 40 K atoms at Shanxi University and of 6 Li atoms at MIT. 



I. INTRODUCTION 

Thanks to the high-controllability of ultracold atoms in 
interatomic interaction, geometry, purity, atomic species 
and lattice constant (of optical lattices), ultracold atomic 
gases have already become one of the most important 
footstones in modern physics, and can be used in sim- 
ulating various strongly correlated many-body models 
in solid state physics Very recently, another con- 

trollability of ultracold atoms is realized experimentally, 
in which the spin degree of freedom of atoms (i.e., the 
atomic internal hyperfine states) is coupled to the or- 
bital degree of freedom (i.e., the momentum) by using 
a two-photon Raman process. This so-called spin-orbit 
(SO) coupling was first created and detected in atomic 
87 Rb Bose-Einstein condensates (BEC) [J, [5|, and then 
produced in atomic Fermi gases of 40 K atoms M and 6 Li 
atoms Q. The realization of SO coupled atomic gases 
enables the simulation of charged particles by using neu- 
tral atoms, which are cleaner and more controllable Q. 
These experiments open an entirely new way to study 
the celebrated effects of SO coupling, such as topological 
insulators, topological superconductors and new exotic 
superfluid phases |9l— flll| . 

Radio-frequency (rf ) spectroscopy has been widely ap- 
plied in many experiments to study fermionic pairing in 
a two-component Fermi gas near Feshbach resonances 
when it crosses from a Bardeen-Cooper-Schrieffcr (BCS) 
superfluid of weakly interacting Cooper pairs over to a 
BEC of tightly bound molecules (l2l . Il3| , and also to in- 
vestigate the properties of polarons in spin-imbalanced 
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Fermi gases This technology allows a precise de- 

termination of the molecular binding e nerg y 
pairing gap in a degenerate Fermi gas [l2| |. In addition, 
the momentum- resolved rf spectroscopy, i.e., the spec- 
troscopy at a specific momentum, is also available (l9| 
and provides important information on low-energy exci- 
tations of novel exotic state of matter in ultracold atoms. 
In this respect, it would be a powerful tool for charac- 
terizing the recently realized spin-orbit coupled atomic 
Fermi gases. Indeed, momentum- resolved rf spectroscopy 
of non-interacting spin-orbit coupled 40 K and 6 Li atoms 
has already been reported [g, [Z| . 

In this work, we aim to qualitatively predict the 
momentum-resolved rf spectroscopy of a resonantly in- 
teracting atomic Fermi gas with equal Rashba and Dres- 
selhaus spin-orbit coupling, a system already realized ex- 
perimentally at Shanxi University and also at MIT by 
using broad Feshbach resonances [5^, [2l|. As is well- 
known, a strongly interacting Fermi gas near Feshbach 
resonances is notoriously difficult to handle theoretically 
[22l. |23| . Therefore, at the qualitative level, we approx- 
imate the strongly interacting Fermi gas as a mixture 
of non-interacting atoms and molecules. All the prop- 
erties of individual molecules are determined from the 
two-particle solution of two-interacting atoms. Our ap- 
proximation may be justified at large temperatures well 
above the superfluid transition temperature T c , where 
molecules are formed below the Feshbach resonance and 
have little correlations among themselves or with atoms. 
This is exactly the situation in current experiments, for 
instance, at Shanxi University @, the typical temper- 
ature of spin-orbit coupled 40 K atoms is now at about 
0.6Ti?, where Tp is the Fermi degenerate temperature. 

We consider separately the momentum-resolved rf- 
spectroscopy of atoms and of molecules. Moreover, ac- 
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cording to real experiments we take into account the 
crucial trapping potential that prevents the atoms and 
molecules from escaping. Within local density approx- 
imation (LDA), the trapped Fermi cloud may be di- 
vided into many locally uniform subsystems. For each 
subsystem, we calculate the rf responses of the atomic 
and molecular components based on the one-particle and 
two-particle solutions of uniform systems and Fermi's 
golden rule. We finally add up all the local contribu- 
tions to determine the total momentum-resolved rf spec- 
troscopy. We note that due to SO coupling, weakly 
bound molecules with anisotropic mass and wavefunction 
may be formed [24U26| . The bound to free rf transition of 
weakly bound molecules at rest has been predicted [27l |. 
In this work, we will consider the rf transition of weakly 
bound molecules at motion with arbitrary center-of-mass 
(COM) momentum. 

The paper is arranged as follows. In the next section, 
we consider the situation of current experiments, and 
introduce the model Hamiltonian responsible for equal 
Rashba and Dresselhaus SO coupling. In addition, we 
introduce the LDA formulism and present the calcula- 
tion for the chemical potential. Then, the single-particle 
and two-particle problems of a local uniform subsystem 
are solved in Sec. Mil The general formulas for the bind- 
ing energy and wavefunction of two-particle bound state 
with non-zero COM momentum are also provided. In 
Sec. IIV1 we derive the momentum-resolved rf transi- 
tion signals for non-interacting atoms and molecules, re- 
spectively, and then obtain the total rf spectroscopy for 
the harmonically trapped ideal gas mixture of fermionic 
atoms and bosonic molecules. Finally, our main results 
are concluded in SecfVl 



The second term in Ho describes the Raman coupling 
between two spin states with strength , kn is the 
wave-number of two Raman laser beams and therefore 
2hkn is the momentum transfer during the two-photon 
Raman process [6|,lZ|. The interaction strength is denoted 
by the bare interaction parameter Uq, which should be 
regularized in terms of the s-wave scattering length a s , 
i.e., 1/Uq = m/ (ATrh 2 a s ) - (l/V)J2 k m/ (h 2 k 2 ). By us- 
ing Fcshbach resonances, the s-wave scattering length a s 
could be arbitrarily tuned, and the system may undergo 
a crossover from a BCS superfluid of weakly interact- 
ing Cooper pairs to a BEC of tightly bound molecules. 
Therefore, the atomic chemical potential fi may decrease 
as a result of decreasing atomic population, when weakly 
bound molecules are formed. 

As argued in the Introduction section, we will treat the 
Fermi cloud as a mixture of non-interacting atoms and 
molecules, although the degree of freedom of molecules 
is not made explicitly in our model Hamiltonian from 
the outset. We will take the molecule as a two-particle 
bound state and determine all its properties by using 
two-particle solution of two interacting atoms. In the 
harmonic trap, we assume that these molecules experi- 
ence the same harmonic trap as atoms, but with double 
mass, i.e., Vm (r) = m(uj 2 x 2 + uj 2 y 2 + <jJ 2 z 2 ). Note that, 
our treatment of ideal mixture of atoms and molecules 
becomes exact for a Feshbach resonance with zero reso- 
nance width. However, for the broad Feshbach resonance 
used experimentally, this treatment is valid at the quali- 
tative level only. 

B. Local density approximation 



II. MODELS 



A. Hamiltonian 



We consider a SO coupled Fermi gas with atomic mass 
m in a harmonic trap V (r) = m(u 2 x 2 +co 2 y 2 + uj 2 z 2 )/2. 
The SO coupling is created by a two-photon Raman pro- 
cess 0, 0|- The many-body Hamiltonian responsible for 
this process may be modeled as H = Ho + "Hint , where 



Ho — 



dr^l (r 



h 2 V 2 
2m 



+ V(t)-h 



dr 



(r) + 
(1) 



is the single-particle Hamiltonian and 

H in t = U J dr*\ (r) *\ (r) * i (r) * t (r) (2) 

is the interaction Hamiltonian describing the contact in- 
teraction between two spin states. Here, 'J/J. (r) is the 
creation field operator for atoms in the spin-state a. 



If the number of atoms is sufficiently large, it is rea- 
sonable to assume that the trapped cloud may be di- 
vided into many locally uniform subsystems with a local 
chemical potential [28l. |29|. Then, within LDA, the ex- 
ternal trap V (r) in the Hamiltonian |lj is absorbed into 
the chemical potential, and we can define effective local 
chemical potential 
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M ( r ) = A* - + "ty + utz 2 ). 



(3) 



Note that the global chemical potential /i, which is de- 
termined by the total number of atoms, can be regarded 
as the local chemical potential in the trap center (r = 0). 

In order to evaluate the chemical potential /i, let us 
consider the ideal gas mixture of fermionic atoms and 
bosonic molecules at finite temperatures. The distribu- 
tions of the atoms and molecules are respectively deter- 
mined by the Fermi-Dirac and Bose-Einstein distribu- 
tions, 



f A (e) 



fu (e) 



1 



=(e-f<)/(faTi) + 1' 



1 



e (c-/-iM)/(fcsT M ) _ I '■ 



(4) 



(5) 
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where ks is Boltzmann's constant, and Ta,m, P and Pm 
are the temperatures, chemical potentials of atoms and 
molecules, respectively. From the thermal and chemical 
equilibrium conditions, we have the following relations, 
Ta = Tm = T and pM — 2p + eb, where sb is the 
binding energy of the molecules relative to the threshold 
of two free atoms. To ease the numerical workload, we 
will approximate it as eb — Eb = fi 2 / (ma 2 ). Here, Eb 
is the binding energy in the absence of SO coupling. 

The populations of the atoms Na (in a single spin com- 
ponent) and noncondensed molecules Nm are given by 

/>oc 

N A = / dtp a (e) I a (e) , (6) 
Jo 



Then from Eq. (|13[) . the chemical potential p can nu- 
merically be evaluated at a specific temperature and in- 
teraction strength Eb- The molecule BEC occurs below 
the critical temperature, i.e., T < T c . In this case, the 
chemical potential of molecules vanishes, i.e., pm = 0, 
which in turn gives p = —Eb/2. 

With the chemical potential p obtained, the properties 
of the trapped systems may be calculated by integrat- 
ing over the whole space, based on the local solutions 
of uniform systems. In the following sections, we will 
firstly solve single-particle and two-particle problems in 
uniform systems, and then investigate the rf spectroscopy 
of trapped gases under the LDA approximation. 



N 



depM (e) hi (e) , 



(7) 



where pa.m (e) are the density of states (DOS) of the 
atoms and molecules, respectively. Qualitatively, for 
atoms we shall use the well-known expression for DOS 
in harmonically trapped systems without SO coupling, 



p A (e) de 



2{hw) 



rde, 



(8) 



where w = (u} x ojyU> z )^ 3 . For molecules, we instead use, 



PM (e) de = 3 <*e- 
2 (rkj) 



(9) 



The factor appearing in the DOS of molecules is due 
to the SO coupling, which induces an effective mass M x = 
7 • (2m) as we will make clear later. Combining Eqs. @ 
and and the populations of the atoms and 

noncondensed molecules are given by 



N A =- 



Ei-i {-za) , 



Nm = V7 f -r— - J Li 3 (z M ) 



(10) 



(11) 



where z A = e^^ kBT) and z M = e ( 2 ^+ E B)/ik B T) are the 
fugacities of the atoms and molecules, respectively, and 
Li n (z) is the polylogarithm function. Therefore, the 
chemical potential p should satisfy the following equa- 
tion, 



2N A 



N m + N$ =N, 



(12) 



where is the population of the condensed molecules, 
and N is the total number of atoms. Note that one 
molecule is counted as two atoms. Above the critical tem- 
perature, i.e., T > T c , there is no condensed molecules 
(N^g = 0), therefore Eqs. yield 

(hcu \ ^ 



III. SINGLE- AND TWO-PARTICLE 
PROBLEMS FOR UNIFORM SYSTEMS 

A. Single-particle solution 

The single-particle problem in a uniform system has 
been discussed in detail in our previous work (27L l30l |. 
Here, for self-containment we summarize briefly the re- 
sults. We focus on the following non-interacting Hamil- 
tonian for the subsystem, 

Ho = E/ dr*t(r)^* a (r) + 

dr h\ (r) e akRX ^ i (r) + ff.c.l , (14) 
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where k 



-iV. Under the transformation 1271 \3£ 



* T (r) = e +tkRX ipt (r) and 4^ (r) = e - lkRX ipi (r), the 
spatial dependence of the Raman coupling term in 
Eq. (fl4|) is removed and yields 



Ho 



dr-0j (r) 
drip] (r) 



h 2 (k 2 + k%) 

2m 

h 2 (k 2 + k 2 R ) 

2m 



Xk x 



\k x 



( r ) + 
^ (r) + 



rfr 



V4 (r) (r) + H.c. 



(15) 



where e x is the unit vector along x direction. For con- 
venience, the effective Zeeman field h = £Ir/2 and the 
SO coupling strength A = H 2 ku/m are introduced. This 
single-particle Hamiltonian (|15D can be easily diagonal- 
ized to yield two eigenvalues |27l. |30T | , 



2m 



h 2 k 2 

2m 



± y / h 2 + X 2 kl 



(16) 



The symbols "±" stands for the two helicity branches, 
and the corresponding single-atom eigenstates in helicity 
basis take the form 



Ik- 
Ik- 



cos 6*k sin 6>k 
— sin 0k cos 0k 



|kt) 

in) 



(17) 
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where 9^ = arctan {^\/h 2 + X 2 k 2 — Afc x J /h . It is ob- 
vious that the lowest single-particle energy occurs at 
k± = \Jk 2 + k 2 = and k x 
h < m\ 2 /h 2 , and is given by, 



y /m 2 X 2 /h i -h 2 /X 2 if 



2 1,2 



Erain — 



h l k 



mX 



2 1,2 



K 2 h 



2 h 2 



2m 



2h 2 2mA 2 



2mA 2 ' 



(18) 



B. Two-body bound state 



The two-body problem with zero COM momentum in 
a uniform system has already been solved in [27} • How- 
ever, real experiments are carried out well above the crit- 
ical temperature [6j, [D] , where two-particle states with a 



non-zero COM momentum may become important. As 
we shall see, at the quantitive level, molecules with a 
nonzero COM momentum can have small contribution 
to the (momentum-resolved) rf spectroscopy. Here, we 
outline briefly the general solution with arbitrary COM 
momenta. 



Under the gauge transformation for the field operators, 
the form of the interatomic interaction @ is invariant, 

Hint = U J drip\ (r) (r) ip± (r) V>t (r). For a two-body 
problem in the presence of SO coupling, the COM mo- 
mentum is a good quantum number, and therefore the 
two-body wavefunction at a given COM momentum q 
can generally be written as (we set the volume V = 1) 

M 



* 2s (q))= i^K (k)c ! 



q/2+kf-q/2-14 



q/2+k|"q/2-kt 



Iptf ( k ) C q 



cLn + ^44 (k) cL„ ,. r 



q/2+kt°q/2-kt 



q/2+k4. q/2-k; 



(19) 



where and c£, are creation field operators of spin-up and spin-down atoms with momentum k, and C = 

EkH^Tl (k)| 2 + IV'J.t (k)| 2 + IV'tt (k)| 2 + IV'-U (k)| 2 ] is the normalization factor of the two-body wavefunction. The 
Schrodinger equation % \^2B (q)) = £b (q) 1^25 (q)) is equivalent to the following set of coupled equations, 



£b (q) - 
£b (q) - 
£b (q) - 
£b (q) - 





ft 2 *; 2 
h 




m 


m 


4m 




fi 2 /fc 2 
h 




m 


m 


4m 


A 2 4 


h 2 k 2 
+ 




\ m 


m 


+ -r~ - 

4m 




ft 2 /fc 2 
+ 




\ m 


m 


+ -r~ - 

4m 



2Afc 2 
2Afc 2 
- Xq a 
Xq, 



i>n ( k ) 
VM-t ( k ) 
V'tt ( k ) 
V^u ( k ) 



+(c/ /2) ]T k , [V'n ( k ') - ( k ')l + ^tt ( k ) + W) 

-(Ua/2) [<p n (k') - V^ t ( k ')l + hiPn ( k ) + W U (kfel) 

(22) 



ft^4. (k) + (k) 
/i^ ti (k) + htl>x (k) 
I 



Here, eb (q) is the energy of the two-body wavefunction 
which is dependent on the COM momentum q. By in- 
troducing 



(23) 



^4kq = £b (q) 



h 2 k\ 

m 



h 2 k 2 h 2 q 2 



rn 



4m 



(24) 



^ (k) = ty n (k) - ^ (k)]/V2 and % (k) = ty n (k) 
-0^t (k)]/\/2, h is easy to show that, 



ip n (k) 
4>U ( k ) 

V>a (k) 



y/2h 



^4kq — A(7 X 

s/2h 
Akq + Xq x 
2Xk x 



i>a (k) 
i>a (k) 



4fr 2 + A 2 9 2 

(x 2 q 2 x , 4/i 2 A kq 



^ ^-^ 2 -X 2 q 2 



(25) 
(26) 

^ (k) ,(27) 
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and 



4A^|A4kq ^(k) = [/ o y> s (k') 



1 - 4/iV (A 2 q - 

From Eq. (|28p . the un- normalized spin-singlet wavefunc- 
tion ip s (k) is given by, 



i> s (k) 



4A 2 fc2/A Uq 



1 - ^ 2 / (^ q " A 2 <z 2 ) 



(29) 



Thus, the two-body binding energy is determined by 
(7oE k ^(k) = l,or 



4tt?i 2 



0. 



(30) 



Here, the bare interaction Uq has been regularized by the 
s-wave scattering length a s . A bound state exists if its 
energy satisfies eb (q) < 2E min , where E min is the lowest 
single- atom energy given by Eq. (|18p . The normalization 
factor for the wavefunction is given by 
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Figure 1: (Color online) Binding energies —eb/Ex of two 
atoms with zero COM momentum q = as functions of the s- 
wave scattering length 1/ (kna s ) for three typical values of the 
Zeeman field. The horizontal dotted lines are the threshold 
energies 2E m in where the bound states disappear. Here, we 
have introduced kn — mX/h 2 and E\ = mX 2 /ti 2 as the units 
of the wavevector and energy, respectively. 



c = ]ThMk)| 2 <h 



4A 2 fc 2 /A 2 



Ah 2 A 2 + A 2 9 2 



A 2 
kq 



A 2 <Z 2 



kq 



1 



4^7 (^q-A 2 ? 2 ,) 



(31) 



The bound state energy can be easily solved from Eq. 
([3H)l and the results are plotted as a function of the s- 
wave scattering length in Fig[T] for zero COM momen- 
tum at three typical values of the Zeeman field (h/E\ = 
0.2,0.5,0.8). Here, we have introduced kn = mX/h 2 and 
E\ = m\ 2 /H 2 as the units of the wavevector and energy, 
respectively. We find that the bound state only exists for 
a s > 0. 

Due to the SO coupling, the effective mass of the two- 
body bound state is affected. In the limit of small COM 
momentum, i.e., q rs 0, the energy of two-body bound 
state would have a well-defined dispersion, 



e B (q) = e B (0) 



2M X 2M X '' 



(32) 



where eb (0) is the two-body binding energy with zero 
COM momentum, and M Xj ±_ are the effective masses in 
x and transverse (y and z) directions, respectively. The 
transverse effective mass M± is not affected by the SO 
coupling, i.e., Mj_ = 2m, since the SO coupling is only 
applied along x direction, therefore 



^kq = A v 



h 2 q 2 ( 2m 



4m V M x 



- 1 



(33) 



where = eb {§) — h 2 k 2 l /m — ti 2 k 2 /'m. Substituting Eq. 
(|33| into Eq. ([30]) . expanding up to the second order of 
q Xl and comparing the corresponding coefficients, we can 
easily obtain 



1 

7 



2m 
M~ x 



= 1 - 



h 2 



h 



where 



h = - 



4h 2 ) 2 +4X 2 k 2 x (A 



Ah 2 



AA\ {Al-Ah 2 -AX 2 k 2 x Y 
\%X 2 h 2 kl 



k 



AUAl 



Ah 2 - AX 2 kl) 



(34) 

(35) 
(36) 



The effective mass ratio 7 = M x j (2m) as a function of 
the interaction strength is shown in Figl2]for three typical 
values of the Zeeman field (h/E\ = 0.2, 0.5, 0.8). We find 
in the BEC limit, the effective mass is independent of 
the values of the Zeeman field, and 7 approaches unity, 
since two atoms are deeply bound and form a tightly 
bound molecule. When the interaction is tuned towards 
the resonance, the two-body bound state becomes looser, 
and the SO coupling will induce a relatively large effective 
mass. However, the large SO coupling could also destroy 
the bound state at a critical s-wave scattering length as 
we have seen in FigJTJ 



IV. RADIO-FREQUENCY SPECTROSCOPY 

In this section, we will firstly present the basic idea of 
the rf transition based on the Fermi's golden rule, and 
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Figure 2: (Color online) The effective mass of weakly bound 
molecules, 7 = M x / (2m), as functions of the interaction 
strength 1/ (kna s ) for three typical values of the Zeeman field. 



then derive the Frank- Condon factor for the rf transition 
of atoms as well as weakly bound molecules. Finally, un- 
der LDA, we theoretically investigate the rf-spectroscopy 
of a harmonically trapped ideal gas mixture of fermionic 
atoms and bosonic molecules, which is observable in cur- 
rent experiments. 



A. Radio-frequency transition and the Fermi's 
golden rule 

The basic idear of the rf transition is simple. An rf 
field is applied to an atomic Fermi gas with two hyperfinc 
states (denoted as |1) = \\) and |2) = \\)), and drives 
one of the hyperfine states (i.e., \D) to an upper state |3) 
with a bare atomic hyperfine energy difference /iu^ due 
to the magnetic field splitting. The Hamiltonian of this 
rf coupling may be written as, 



Vrf =V dr 



v4 ( r ) *i ( r 



H.c. 



(37) 



where ^3 (r) creates an atom in the third state at position 
r, and Vq is the strength of the rf drive and is related to 
the Rabi frequency lur by Vq — Hujr/2. Using the field 
operator after the gauge transformation, the rf coupling 



Vrf = V I dr 



i>l (r)^(r)+#.c. 



(38) 



(39) 



for later convenience. Note that, because of the transfor- 
mation, effectively there is a momentum loss of kjie x for 
the transferred atoms. 



or in the momentum space, 

Vrf = V J2(<l-k R e„, 3 CH + H.C 



According to the Fermi's golden rule, the general 
strength of the rf transition is proportional to the Franck- 
Condon factor [31], 



r(n) = \(^ f \v rf \^)\ 2 5 n 



E f - Ei 



(40) 



where \tpij) denote the initial and final states with cor- 
responding energy Eij , respectively, and f2 is the fre- 
quency of the rf field. The delta function keeps the energy 
conserved during the transition. 



B. The single- particle rf transition 

In the single-particle picture, the eigenstates of a SO 
coupled atom are the helicity states |k+) (|k— )), other 
than the original spin states |k f) (|k \)). If an atom is 
initially prepared in the state |k— ) with energy E k _, the 
rf photon will transfer this atom to a third empty state 
1 3). In order to obtain the final state of the rf transition, 
it's useful to calculate V r /|k— ). Using Eqs. (jTTJ) and 
(]39| . we can easily obtain 



V 



rf 



V -cos6 k ct_ kRt 



\vac) 



(41) 



This means after the rf transition, there is a probability 
of I cos 8k I that transfers the atom to the final state |3) 
with a momentum k — k^e x . Similarly, the final state 
should be 



Vrf |k+) = U c 



-sinf/kc^^ 



(42) 



if the atom occupies the state |k+) before the rf tran- 
sition. It also means the probability that transfers an 
atom to the final state |3) with a momentum k — kne x 
is |sin#k| 2 - Generally, the atoms are initially prepared in 
the mixture states of |k— ) and |k+) with the probabil- 
ity / (E k ± — m) i respectively, where / (i?k± — A*) is the 
Fermi-Dirac distribution, and /x is the chemical poten- 
tial. According to the Fermi's golden rule, the total rf 
transition strength therefore is given by (Vq = 1) 

r A (fi) = ]T sin 2 e k f (E k+ - fi) s (n - ^±J + 



cos 2 e k f(E ]i --fi)s(n 



Ac 



,(43) 



where 



Ae± = W^H E k ± 

2m 



= -Xk x T Vh 2 + A 2 fc2 



(44) 



is the probable energy difference between the final and 
initial states. Without any confusion, we have ignored 
the bare hyperfine splitting, i.e., uj^i = 0. 
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C. Two-body bound to free radio-frequency 
transition 

The picture of the two-body bound to free rf transition 
is that the rf photon breaks a weakly bound molecule and 
transfers one of the two atoms to the third state |3). The 
final state is determined by [13] 

V rf \*2B (q)> =V J2 4- fcB e„ 3 Ck'| \*2B (q)) • (45) 



Substituting the wavefunction of the two-body bound 
state (fT§|) into Eq. (|4"5"]l . and after some straightforward 
algebra, we shall easily obtain 



-'q/2-k-fc. R e x ,3 / 



Vrf\V 2 B (q)) = -yf^oE< 

(k) + <p a (k)] Cq/2+kt + Vty u ( k ) </ 2+ h} 



(46) 



Using Eq. (|17p , the final state can be written in terms of 
the helicity basis as 

V r /|* 2B (q)> = -yiVoE c i/a-k-W x 

k 

( S q/2+k+C q/2+k+ - S q /2+k-C q /2+k-) \ vac ) t (47) 



where 

Sq/2+k+ = [V^s (k) + V^a (k)] COS 6> q/2+k + 

%/2^ u (k)sin0 q/2+k , (48) 

•Sq/2+k- = fya (k) + Va (k)] sin 6» q/2+k - 

V^Vu (k) cos # q/2+k . (49) 

Eq. (|4"T)) could be understood as follows. During the rf 
transition, a molecule with COM momentum q is bro- 
ken by absorbing the rf photon, and then the spin-down 
atom is transferred to the third state. Finally, there are 
two possibilities. We may have two atoms in the third 
state and upper helicity state, respectively, with a prob- 

l 1 2 

ability |s q / 2+k+ | /C Also, we may have a probability 

l 1 2 

°f s q/2+k- /C for having two atoms in the third state 
and lower helicity state, respectively. Taking into ac- 
count these two final states, we should have the following 
transfer strength or the Franck-Condon factor, according 
to the Fermi's golden rule (Vq = 1), 





V 




«q/2+k+ ^ 










K/2+k- ^ 6 


h I 



, (50) 



where 
Ae ± = 



h 2 (q/2 k k R e x ) 2 
2m 



■£ q /2+k,± ~e B (q) (51) 



is the probable energy difference between the final and 
initial states. Here, we have already ignored the bare 
hyperfine splitting, i.e., W31 = as well. Eqs. (|50|) and 
(f5Tj) recover the results of [27( when q = . 



Rf spectroscopy for a harmonically trapped gas 
mixture 



D. 



In this subsection, we are ready to investigate theoret- 
ically the rf response signals of the harmonically trapped 
ideal gas mixture of the atoms and molecules, as well as 
the total rf spectroscopy. 



1. The signals of the atomic component 

The rf transition strength for atoms in a uniform sys- 
tem is determined by Eq. (|4"3")h In the harmonic trap, the 
rf transition strength of atoms at the position r should 
be (O, r) , in which we shall use the local chemical 
potential, 



A* (r) = p - j(u 2 x x 2 + UyV 2 + u 2 z 2 ), 



(52) 



where /i is the chemical potential in the trap center de- 
termined by Eq. ([1U1) . Thus, using LDA the total rf tran- 
sition strength of the atomic component in a harmonic 
trap is T A (ft) = / drT A (ft, r). 

In the actual calculation, the summation ^ k 
in r^4 (ft, r) may be replaced by the integral 
(2tt)" 3 / dk x J d 2 k ± . After some tedious but straightfor- 
ward derivation, the total atomic rf transition strength 
has the following form, 



ExT A (») 
N 



2V2vr 



Ep Tp 



dk±Li 3/2 



-z A e 



3/2 h2 
ft 2 "* 

E k /(k B T) 



where 



E k 



h 2 (k 2 ± 



4) 



2m 



2mA 2 



! -ft 2 
2ft 



,(53) 



ft 2 



2 ft 



(54) 

and Ep = hu (3iV) 1/3 and Tp = E F /k B are the Fermi 
energy and Fermi temperature, respectively. At a given 
temperature T and interaction strength Eb, the atomic 
chemical potential fi is obtained from Eq. (|13|) . Then, 
with a fixed total atom number N (which in turn gives 
the Fermi energy Ep), the rf transition strength of the 
atomic component can numerically be calculated from 
Eq. (|53p. The results are presented in Fig. [3] for dif- 
ferent values of Zeeman field at the given parameters, 
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Figure 3: (Color online) The rf spectroscopy of the atomic 
component in a harmonic trap at different values of the Zee- 
man field, i.e., h/E\ = 0.6,0.8,1.0. Here, we have taken 
Eb/Ex = 2, E F /Ex = 1, and T/T F = 1. 



E B /E X = 2, E F /E X = 1 and T/T F = 1. We find that 
there is a bimodal structure in the rf spectroscopy of the 
atomic component. The stronger rf signals at the posi- 
tive frequency Q correspond to the transition out of the 
lower helicity state, while those with weaker strength at 
the negative frequency correspond to the transition out 
of the upper helicity state. On the other hand, there is 
a red-shift in the peak positions as the strength of the 
Zeeman field decreases. As we may anticipate, the peak 
transition frequency Q will reach the bare hyperfine split- 
ting W34. = in the absence of SO coupling. 



In addition to the total strength of the atomic rf tran- 
sition, we can also calculate the transition strength of 
atoms at a specific momentum k x . This momentum- 
resolved rf spectroscopy has already been investigated 
at Shared University [a] and at MIT [7]. According to 
Eq. ([43]) . we can carry out the integral over kj^ and re- 
tain that of k x . Therefore, the momentum-resolved rf 
transition strength is determined by 



J 



k R E x T A (Q,k x ) = _30r (Ex _ T_y /2 { | _ 



Xk r 



2V2 \E F T f 

VtF /Ex jt \ 3/2 / Afc, 

2V2 \e f 't f ) 



\/h 2 + X 2 k 2 



^h 2 + \ 2 k 2 x 



5 fi- 



8 n 



h 



roc 

J dklLi m ^z A e- E ^' kBT 



dk 2 Li 



-E k _/k B T^ 



(55) 



Corresponding to experiments, the delta function in Eq. 
([S3)) can be replaced by the Lorentz-line shape as, 



where a is the line width. 

Fig. [4] reports the momentum- resolved rf spectroscopy 
of the atomic component for different values of the Zee- 
man field and interaction strength. Here, we define 
K x = k x — ktf to make sure the spectra is a symmet- 
ric function of K x in the absence of SO coulping (30j . 
As we anticipate, there are two branches in which one is 
relatively weaker than the other one. Since the energy of 
the upper helicity state |k+) is larger than that of the 
lower helicity state |k— ) , the atoms will occupy the lower 
helicity state first and then the upper one before the rf 
transition. Thus, the initial atomic population of the 
lower helicity state should be larger. Therefore, similar 
to the integrated rf spectroscopy, the brighter branch of 
the momentum-resolved rf spectroscopy corresponds to 
the rf transition out of the state |k— ), while the weaker 
one corresponds to the rf transition out of the state |k+) . 
Therefore, the contribution from two initial states could 
be well identified experimentally. We can also see that 



the gap between two branches becomes larger when the 
strength of the Zeeman field increases. This is because 
the transition frequency deviates from the bare hyperfine 
splitting more obviously with a stronger SO coupling. 

All these features, discussed in the above, are in good 
agreement with the experimental observation for a non- 
interacting spin-orbit coupled atomic Fermi gas: see, for 
example, Fig. 4 in Ref. Q and Fig. 2 in Ref. 0. 



2. The rf spectroscopy of the molecular component 

With the two-body binding energy solved from Eq. 
([50)) at a given COM momentum q, the rf transition 
strength of a single molecule can be calculated directly 
from Eq. ([STJ)) . We plot the rf response of a single 
molecule with zero COM momentum, as well as non- 
zero COM momentum (i.e., q x /kn = 0.5) in Fig. [5] 
Here, the interatomic interaction strength, denoted by 
Eb/E\ = [h 2 / (m\a s )] , is set to be unity. Since the 
SO coupling is only along the x direction, we can set the 
transverse COM momentum qj^ = 0, which is physically 
irrelevant and arbitrary. From Fig. [5) we find that there 
is a slight red-shift in the peak position with a non-zero 
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Figure 4: (Color online) The momentum-resolved rf spec- 
troscopy of the harmonically trapped atomic component for 
different values of the Zeeman field and interaction strength. 
Here, we have taken E F /E X = 2, T/T F = I and a = O.f. In 
addition, we have also defined K x = k x — kn (see text), which 
moves the whole spectra to the right by an amount ka. 




Figure 5: (Color online) The rf spectroscopy of a single 
molecule for different values of the Zeeman field at given COM 
momenta (\/kn = and q/kn = 0.5. Here, we have chosen 
an interaction strength Eb/E\ = 1. 
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Figure 6: (Color online) The momentum-resolved rf spec- 
troscopy of a single molecule for different values of the Zeeman 
field at given COM momenta q^/fcij = and q/fcji = 0.5. 
Here, we have chosen the interaction strength Eb/E\ — I. 



COM momentum. This is because the energy gap be- 
tween the final and initial states becomes smaller at non- 
zero COM momentum, and consequently a little higher 
transition peak appears. 

The momentum-resolved rf spectroscopy of a single 
molecule can be obtained from Eq. (|50[) by integrating 
over kj_, and is given by 



T M (fi, k x ) 



8ir 2 hC 



K/2+k+| 2 e>(fcL+) + 
K/2+k-f e> (k 2 ± _) 



where O (x) is the step function and 



m 



kl ± = -[Q + e B (q)]- 



k 2 x + (k x + k R y 



^ + (fc f + g - )2 ±^^ + A^ ( q ^ + k x 
4 r V V 2 



(57) 



(58) 



The momentum-resolved rf spectroscopy of a single 
molecule with zero and non-zero COM momentum 
(ix/kR = 0.5) are presented in Fig. [SJ In this plot, we 
have defined K x = —k x — kji, as we require that in the 
absence of SO coupling the spectra is an even function 
of K x [27} • We find that the contribution from two final 
states are well separated in different frequency domain 
and therefore should be easily observed experimentally. 



3. The total rf spectroscopy of the gas mixture of atoms 
and molecules 



By increasing the attractive interaction strength of a 
two-component Fermi gas at finite temperatures, the sys- 
tem will evolve into a mixture of fermionic atoms and 
bosonic molecules. We have so far discussed the in- 
tegrated and momentum-resolved rf responses for both 
atoms and molecules. The total rf spectroscopy of a 
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Figure 7: (Color online) The total rf spectroscopy of a har- 
monically trapped ideal gas mixture of fermionic atoms and 
bosonic molecules at different Zeeman field. Here, we take 
E B I Ex = l,Ejp/Ex = \ and T/T F = 0.2 . 



Figure 8: (Color online) The total rf spectroscopy of a har- 
monically trapped ideal gas mixture of fermionic atoms and 
bosonic molecules as a function of temperature at a fixed Zee- 
man field strength h/E\ = 0.5. Here, we use Eb/E\ = 1 and 
E F /E X = 1. 



harmonically trapped ideal gas mixture of atoms and 
molecules can now be easily obtained. 

For the atomic component, the rf transition strength 
is already given by Eq. f|53[) . As to the molecules, the 
population Nm can be calculated from Eq. (fTTj) at given 
temperature T/Tp and interaction strength Eb/E\. In 
real experiments, the temperature is low so that the dis- 
tribution of the molecular COM momentum q is very 
narrow around q = 0. As already seen from Figs. 5 and 
6, the rf spectroscopy of molecules depends weakly on 
q. Thus, as a good qualitative approximation, we may 
focus on the rf transition of molecules with zero COM 
momentum only. Then, the total rf transition strength 
of molecular component is given by NmTm (^)- hi the 
end, we obtain the total rf spectroscopy, 



V. CONCLUSIONS 

In conclusion, we have theoretically investigated the 
momentum-resolved radio-frequency spectroscopy of a 
harmonically trapped ideal gas mixture of fermionic 
atoms and bosonic molecules with spin-orbit coupling. 
This is a qualitative theoretical model of spin-orbit cou- 
pled atomic Fermi gases near broad Feshbach resonances. 

In particular, for local uniform systems, general for- 
mulas for the wavefunction and binding energy of weakly 
bound bosonic molecules with non-zero center-of-mass 
momentum have been derived. The influence of this non- 
zero center-of-mass momentum on the radio-frequency 
spectroscopy of molecules has been discussed in detail. 



r t (ft) = r A (n) + n m t m (fi) . (59) 

Fig. [7] shows the total rf spectroscopy of the atom- 
molecule mixture at different Zeeman field and at fixed 
temperature T/Tp = 0.2 and interaction strength 
Eb/E\ = 1. It is obvious that there is a double-peak 
structure in the spectroscopy. The peak at higher fre- 
quency is responsible for weakly bound molecules, since 
more energy is required for pair breaking. With increas- 
ing the Zeeman field h, the energy difference between the 
final and initial states of atoms during the rf transition 
becomes larger, thus a decrease in the atomic peak is ob- 
served. On the other hand, since the energy difference 
between the molecular final and initial states becomes 
smaller as the Zeeman field increases, the rf signal of 
molecules grows, as we may anticipate. 

In order to show the temperature dependence, we also 
calculate the total rf response of the ideal gas mixture 
at different temperatures and at a fixed Zeeman field 
h/E\ — 0.5, as reported in Fig. [5] As the temperature 
decreases, the population of molecules becomes larger. 
Thus, a reduction of the atomic signal is observed. The 
molecular response should dominate when temperature 
decreases to zero. 

Finally, we present the total momentum-resolved rf 
spectroscopy in Fig. |H1 with parameters correpsonding 
to Fig. [7J (see the left panel) and Fig. [5] (right panel). 
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Figure 9: (Color online) The total momentum-resolved rf 
spectroscopy of a harmonically trapped ideal gas mixture of 
fermionic atoms and bosonic molecules. The left panel shows 
the spectra at different Zeeman field strength at a fixed tem- 
perature T/Tf = 0.2, while the right panel shows the tem- 
perature dependence of spectroscopy at a fixed Zeeman field 
hf E x = 0.5. Here, we take E B /E X = 1 and E F /E X = 1. 
Note the different color scale in the left and right panels. 
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Based on few-body solutions in uniform systems, the 
radio-frequency responses of atomic and molecular com- 
ponents in harmonic traps have been calculated within lo- 
cal density approximation, respectively, by using Fermi's 
golden rule. Our results provide a qualitative picture 
of momentum-resolved radio- frequency spectroscopy of a 
strongly interacting spin-orbit coupled Fermi gas in har- 
monic traps, which could be directly observed in future 
experiments at Shanxi University [1,[20| and MIT^I^. 

In real experiments, where strongly interacting Fermi 
gases are created in the vicinity of Feshbach resonances, 
the atoms and weakly bound molecules must strongly 
interact with each other. An in-depth analysis of radio- 
frequency spectroscopy would rely on much more com- 
plicated many-body theories [32||, beyond the simple few- 



body picture presented in this work. 
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